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LEXINGTON 


MASSACHUSETTS 


ABSTRACT 


Cramer-Rao  lower  bounds  on  the  degradation  of  measurement 
precision  for  closely  spaced  optical  targets  are  found.  Degra- 
dation is  more  severe  than  previously  thought*  especially,  for 
amall  target  separations.  The  bound  obtained  by  a previous  analysis, 
i»,  in  fact,  theoretically  identical  to  the  Cramer-Rao  bound  pre- 
sented here,  Tho  disagreement  in  results  is  caused  by  numerical 
problems  in  a computer  program  used  in  the  previous  analysis. 
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I.  INTRODUCTION 


During  the  past  several  years  considerable  attention  has  been 


focused  upon  the  accuracy  of  parameter  estimates  for  optically  un- 


resolved objects  [1]  - [41.  A review  of  the  literature  shows 


that  predictions  for  the  variances  of  the  parameter  estimates  are 


obtained  via  different  analytical  aporoaches  using  different 


assumptions  on  the  noise  environments.  In  |1),  a particular 


"resolution  scale"  is  oroposed  in  order  to  relate  the  resolution 


capability  of  a sensor  with  its  measurement  precision  degradation 


factor  due  to  multi-target  interference.  Similar  studies  have 


been  made  for  radar  system  applications,  see  for  example  [5)  and 
[6  1.  It  is,  moreover*,  *uOintad  ftift*li\  "|  5 I that*a‘ como’lete  analysis 


of  the  resolution  problem  requires  extensive  simulation  and  that 


the  resulting  error  probabilities  may  depend  heavily  upon  the  choice 


of  the  detection  as  well  as  the  estimation  Bchemo.  The  criterion 


for  resolving  two  closely  spaced  radar  targets  appears  to  be  some- 


what arbitrary.  In  16 1,  for  example,  it  is  defined  as  twice  the 


value  of  the  angular  accuracy  of  the  weaker  target.  Similar  measures 


have  also  been  used  in  the  optical  community.  It  is  outside  the 


scope  of  this  report  to  discuss  the  merits  of  various  resolution 


criteria  in  a more  extensive  manner.  It  is,  however,  important 


to  noint  out  that  using  parameter  estimation  performance  degradation 


results  from  111  - [4 1,  or  other  similar  works,  to  predict  the 


resolution  capability  of  some  other  specific  sensor  svBtem  without 


extensive  simulation  can  bo  misleading.  The  casual  use  of  numerical 


r- 


results  from  these  works  could  lead  to  false  conclusions  about 
the  sensor  capability. 

In  l 1 1 - [4],  theoretical  lower  bounds  on  the  variances  of 
various  parameters  associated  with  the  target  intensities  and  angular 
positions  have  been  derived.  The  bounds  presented  in  [2]  - [4] 
are  Cramer-Rao  lower  bounds.  Different  pulse  shapes  and  sensor 
noise  models  are  considered  and  compared.  The  results  presented 
in  (11  are  obtained  using  a different  error  analysis  technique 
and,  furthermore*,  the  parameters  to  be  estimated  are  different 
from  those  considered  in  |2|  - [4].  Attempts  have  been  made  by 
Miller  [2]  to  compare  results  obtained  from  the bc  two  different 
techniques  and  problem  formulations.  Unfortunately,  Miller  could 


not  achuive  agreement  between  the  numerical  results  nor  explain 
the  disagreements.  Miller's  predicted  degradation  Df  the  parameter 
estimation  performance  was  much  more  Bovere  than  that  presented 
by  Pried  [1],  particularly  for  target  separations  much  less  than 
the  optical  diffraction  limit.  Surnrisingly  good  'Stimation  per- 
formance was  predicted  in  f 1 ] for  small  target  separations.  These 
results  have  been  used  on  several  occasions  for  oredicting  re- 
solution capability  of  certain  optical  system.  For  this  reason 
we  felt  it  important  to  establish  the  cause  of  the  discrepancies 
reported  in  [2], 

As  discussed  in  (21,  although  the  problem  formulations  are 
quite  different  in  III  and  |2|,  they  should  produce  results  which 
are  comparable.  The  two  analysis  differ  in  the  choice  of  ootical 
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pulse  shape  but  we  do  not  believe  that  this  alone  could  be  the 
cause  of  the  drastic  discrepancies  observed  in  [?.].  The  purpose 
of  this  report  is  to  establish  the  cause  of  the  above-mentioned  dis- 
agreement. Our  analysis  has  Bhown  that  the  two  approaches  although 
superfically  different  are  in  fact  identical  analytically.  We  can 
demonstrate  that  the  differences  in  results  are  caused  by  numerical 
problems  which  exist  in  one  of  the  computer  programs  used  in  [1], 

In  the  next  section,  the  Cramer-Rao  bounds  are  derived  for  the 
problem  formulated  by  Fried  [1],  A detailed  derivation  which 
established  that  the  equations  in  [1]  are  in  fact  Cramer-Rao  boundB 
may  be  found  in  the  Appendix.  In  Section  III,  the  cause  of  the 


numerical  problems  in  [1] 

% « • jp  m • <«  • k i » ji  v i 


is  discussed  and  illustrated  by  a particular 

. « * • ft  it  * 4 « • 4 s • • * i • » 4 41  . a # . I 4b,  , 


example.  Results  of  some  typical  cases  considered  in  [1]  are  presented 
and  compared.  A summary  and  conclusions  are  given  in  Section  IV. 
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II. 


ANALYSIS 


Let  us  assume  that  p(t)  is  the  output  of  an  optical  sensor  to 
a unit  strength  point  source  such  that 


c 

/ 


p“  (t)dt-l. 


(1) 


The  problem  treated  here  concerns  the  measurement  of  a pair  of 
point  sources  with  relative  strength  (l+*iA)  s (l-*jA)  , separation  T 
and  the  location  of  the  midpoint  between  them  at  tQ.  The  output 
of  the  sensor  becomes 


s(t)  « a ( 1+tjA)  p (t-t0+*}T)  + a (1»*jA)d  (t-t0~>jT)  (2) 

where  a is  the  average  strength  of  the  point  sources  which  does 

•f*  *f 

not  appear  in  the  problem  formulation  given  in  tl).  We  wish 
to  determine  these  parameters  from  a noisy  measurement  taken  at  the 
output  of  the  sensor 


y ( t ) “ s(t)  + n(t) 


(3) 


where  n(t)  is  white  Gaussian  noise  with  two-sided 

fPo'r  easier  comparison  with  results  presented  in  11]  » we  adopt  the 
notation  uned  in  [11. 

ttFor  the  analysis  presented  in  ill,  the  value  of  a is  not  important. 
The  reason  for  having  it  in  Eg.  (2)  will  become  clear  in  later  dis- 
cussions . 
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V-  . 


power  spectrum  density  NQ/2. 

Let  p < t ) be  the  autocorrelation  function  of  p(t),  i.e. 


p (t ) p (t-t) dt 


and  p(0)=l.  The  Cramer-Rao  lower  bound  on  the  variance  of  any 
unbiased  estimator  of  the  unknown  parameters:  relative  strength, 

R(*  (1“*jA)  / (1+JjA)),  separation,  T,  midpoint,  tQ,  can  be  obtained 
by  inverting  the  Fisher  information  matrix,  F,  whoee  (i,j)th  element 

4-.  [7] 


9 1 n A 3 In  A 


where  E{  } denotes  statistical  expectation  and  a ^ denotes  the  i 
unknown  parameter,  namely,  «1»R,  ot^T,  a3»to  and  a4»a,  and  where 
In  A ia  the  log  likelihood  ratio  [7,  p.  274] 

In  ^ “ fT  Jtf  y(t)s(t)dt  - J s2(t)dt|, 


Substituting  (6)  into  (5),  we  have 


, « 2 j r h s (t)\  /3,  s ( t)\ 

lj  *o\l  {>  H)\,  *,) 


where 


3 B — = ‘jafl+SA)2  (p(t-t  “JjT)  - p(t.~trt+JjT)  ] 

a r ° ° 


3 a^---  = Hu  f(l+*5A)p{t-t  +>JT)  - (l-JjA)o(t-t  -*T)1 

5 T ^ ° ° J 

- -a^l+JjA)p(t-t0+»|T)  + (1-»«A)0  (t-t0-*sT)] 


„ (l+>jA)p(t-t  +ST)  + (l-«jA)p(t-t  -VD 
3 a 


and  p(t)  is  the  Is  order  dervative  of  p(t)  with  respect  to  t.  Then 
t-he  elementH  of  the  Fisher  information  matrix  (Eq.  (7))  arei 


F11  - a2  (l+«j A)  4 fl-n(T)  1/NQ 


F12  « SAa2(l+SA)2  p(T)/N0 


F13  - -2a2{l+SA)2  n(T)/No 


F,A  » -Aa  ( 1+*jA)2  [l-p(T)]/Nr 


F22  - a2  [ ( 1-%A2)V<T)  - (1+H A2)  p‘(0)  1/Nq 


F23  - 2Aa2p'(0)/NQ 


• i •<  •••  ■ | « « • «*  « I « • <|  >f  '•>  « 


F 24  = 2a(l-^A^)p(T)/N0 


F 3 3 * -4a2[  (1-*iA2)p‘(T)  + (1+*  A2 ) p’(  0)  ] /NQ 


F34  " 0 


F44  - 4 [ ( l-t-'s  A2  ) + (I-IkA2)  p (T)  1/N0 


where  p and  p*  are  the  1st  and  2n<^  order  derivatives  of  p,  res- 
pectively. 

In  principle,  one  can  invert  this  4:<4  symmetric  matrix,  F, 

d.' VectLy  Vo  *^ind*  tfie*  iower  bounds  on*  the**varia'nce*  of*  the  'parameter' 

estimates  for  oarticular  values  of  R,  T,  t and  a.  if  a is  not  a 

parameter  of  interest,  wo  can  also  invert  an  equivalent  3x3 

matrix,  F,  associated  with  R , T and  t . T can  be  obtained 

o 

from  F by  applying  the  matrix  inversion  formula  for  a partitioned 
matrix  as  below: 


F ' 'V'V  F40_1  P2 


S L 

where  F^  is  the  upper  left  3x3  matrix  of  F,  F2  is  the  1 three 
t h 

elemonts  of  the  4 row  of  F.  Therefore  the  elements  of  F can  be 
wri tton: 


Fu  - 4a2(l+HA)4(l-p2(T)  )/Ff 


(11a) 
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(lib) 


F12  =»  4Aa2  (1+«*A)  2p  (T)/Pq 

F13  = -2a2(l+:SA)2p(T)/N0  (Ho) 

F22  * a2[(l-(,A2)*p*(T)-(l+^A2)'p(0)l/No  -4  {a  (1-*A2)  p (?)  }?/PQ  (Ud) 

r23  * 2Aa2n‘(0)/No  (lie) 

F33  * -4a2(  (l-i<A2)*p(T)+(l+%A2)V(n)  )/No  (llf) 


and  Fg 


dig) 


The  lower  bounds  on  the  variances  of  the  estimates  for  R,  T, 
and  l ares 


CRB 

(R)  = (F  )n, 

-1 

(12a) 

CRB 

(T)  = (F  )22, 

-1 

(12b) 

CRB 

(tQ)  = if  )33. 

(12c) 

The  relative  measurement  precision  defined  in  [1]  can  be  written 
as  the  followings^ 


FA(T) 


A 

rs 


*5 

CRB(R)  J 


(13a) 


TrHiT  relative  measurement  precision  is  the  reciprocal  of  the 
degradation  factor  presented  in  [21- [4).  The  estimation  performance 
degrades  as  the  value  of  the  relative  measurement  precision  decreases. 
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Ft<T) 


A 


CRB(T)Ta 
CRB (T) 


(13b) 


Ft(T) 


CR»(to)T_ 

CRB(t0) 


(13c) 


where  CRB(R)Taoo,  CRB(T)Taa)  and  CRB (tQ) Tbi(B  can  be  obtained  by  in- 
verting F Cor  T®*.  It  is  obvious  that  p (°°)  *M°°)  “p^1”)  then  we  have 


CRB(R)Toi<t,  = NoIl+l(A2l//a2tl+*?A]4, 
CRB(T)T_m  » -N0(l+»<A2jX2|l-><A2lV('l)  , 


(14a) 

(14b) 


and 


CRB  ( t 


o'  T"-1 


» IjCRB(T) 


(14c) 


From  Eqe.  (lla)-(lin  and  (14a)- (14c),  it  can  be  seen  easily  that 
the  degradation  factors  F^f  and  are  in  fact  independent  of 

the  value  of  a.  This  independence  has  also  been  observed  in  (2J- 
1 4 ) for  the  additive  white  Gaussian  sensor  noise  model.  It  should 
be  noted  that  the  actual  estimation  performances  of  the  estimators 
of  R,  T and  t are  not  independent  of  target  intensity  a,  this  can 
be  seen  easily  from  the  expression  of  the  CRB's. 

The  autocorrelation  function  p ( t ) and  its  derivatives  do  not 


appear  specifically  in  the  expressions  for  F.,  F_  and  F.  in  [1] 


irrraB  raars  m 


It  is  not  difficult  to  show  that  they  are  actually  related  to  the 
functions  defined  in  [ 1 1 , namely* 


P (T)  = G 0 (T ) , 


P(T)  » ~2li  Hx  (T>  , 


V(T)  - - (2n)2Gi(t) 


With  the  help  of  Eqs.  (15) -(17),  one  should  be  ablo  to  show  that 
the  degradation  factors  given  in  Eqs.  (13a)-  (13c1)  are  in  fact 
identical  to  those  presented  in  (11.  The  proof  will  be  given  in 
the  Appendix. 


III.  NUMERICAL  RESULTS 


Wo  find  that  the  numerical  values  of  F^,  t'T  and  F^.  calculated 
via  Eqa.  (13a) -(13c)  using  values  of  Gg,  G2  and  provided  in 
Table  1 of  [1)  are  identical  to  values  of  F , FT  and  Ffc  in  Table 
3 of  [1|,  This  proves  that  the  analytical  resultB  of  111  and  the 
Cramer-Rao  bounds  derived  in  the  previous  section  are  identical. 
However,  this  docs  not  explain  the  reason  for  discrepancies  between 
the  results  presented  in  [1]  and  121.  After  we  employ  the  auto- 
correlation  function  of  the  pulse  shape  used  in  [21,  that  is: 


P (T)  « 6 (aT-sinaT)/a3T3 


IThls  Ts  a pulse  shaoe  suggested  by  Fried  [1,  Eq.  (72)]  as  a good 
approximation  to  the  pulse  shape  considered  in  [1], 


• 7’.  . ; • " ■ ..  ■■  . . ..  1 •‘lv  . : i •TTOwmjuiKuiraiHH^fwwdTr-t.c.in 


t v I!- 


with  a«1.6iT,  we  are  able  to  regenerate  results  presented  in  [2] 
by  using  the  computer  orogram  given  in  !1]  - a program  which  cal- 
culates Fa,  Ft  and  Ffc  from  a given  set  of  values  of  GQ,  Gj  and  H^, 
It  is  quite  obvious  that  the  discrepancies  discussed  in  [2]  are 
merely  numerical  problems  which  exist  in  the  program  for  the  cal- 
culation of  Gg,  Gp  and  H^.  It  is  found  that  in  that  program  of 
reference  [l]1  the  following  formulas  were  useds 


sin (T+AT)  *»  sinT  cosAT  + cosT  sinAT 
cos (T+AT)  ■ cosT  cosAT  - sinT  sinAT 


(19a) 

(19b) 


for  numerical  integration  (Simpson's  rule)  of  GQ,  Gj  and 
Interpolation  of  sine  and  cosine  functions  via  Eqs.  (19a),  (19b) 
may  often  have  numerical  inaccuracy,  especially,  when  these 
functions  are  evaluated  at.  multiples  of  it*  We  believe  this  to 
be  the  cause  of  numerical  problem  observed  in  [1],  when  the  target 
separation  is  2.44  X/D.  In  order  to  show  the  effects  of  such 
interpolation  errors  on  results  in  [1],  we  selected  a pulse  shape 
such  that  the  analytical  expressions  for  GQ,  G ^ and  Hi  can  be  obtained. 
The  pulse  shape  having  autocorrelation  function  given  in  (10)  is 
used.  The  numerical  errors  for  0^,  G2  and  Hj.  are  shown  in  Figure 
1 for  this  particular  pulse  shape.  Although  the  percentage  errors 
on  the  values  of  Gn,  G2  and  do  not  appear  significant,  they 
are  large  enough  to  create  severe  discrepancies  in  the  bound 


*Th'e  actual  fault  In  Lhe  calculation  is  not  these  formulas,  but 
the  programming  error  that  updates  the  value  of  sinT  by  sin (T+AT) 
before  calculating  (lib) . 


11 


r 


.2  .4  ,6  ,8  1. 

Target  Separation  In  Units  of  X/D 

<5  value  of  numerical  errors  for  , G ^ and 
er  program  UBed  in  [1]  for  pulae  shape 
elation  function  given  in  Bq.(18)« 


computations.  It  is  easiest  to  see  the  difference  in  the  value 
of  the  determinant  of  the  covariance  matrix 

V A cov(St0,  ST,  <5d]  (20) 

computed  via  interploation  and  exact  values  of  G0,  G2  and 

The  square  roots  of  the  ratioB  of  these  determinants  are  plotted 

against  target  separation  in  Figure  2.  It  indicates  that  there 

are  severe  discrepancies  in  the  values  of  these  determinant!*, 
especially  when  the  target  seoaration  is  small. 

In  the  rest  of  this  report  only  the  Dulse  shape  presented  in 
[ 1 | will  be  considered.  Figure  3 is  a comparison  of  results  ob- 
tained by  the  method  of  this  report  to  results  presented  in  [1] 
and  [2],  As  expected,  results  obtained  from  the  current  analysis, 
or  from  the  numerical  integration  with  exact  values  of  sines  and 
cosines,  are  much  closer  to  those  obtained  via  the  approximate 
pulse  shape  presented  in  [2],  It  is  interesting  to  notice  that 
the  difference  between  current  results  and  results  presented  in 
[1]  is  proportional  to  the  ratio  of  the  square  root  of  the 
determinants  illustrated  in  Figure  2. 

Results  of  F^,  Ft  and  F for  a detector  width  equal  to 
2.0  X/D  from  the  current  analysis  are  shown  in  Figs.  4,  5 and  6, 
respectively.  Values  of  F^,  FT  and  presented  in  [1]  are 
plotted  for  comparison.  We  notice  that  a similar  trend  of  dis- 
agreement is  repeated  here.  For  separation  much  Iobs  than  X/D, 
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Relative  Measurement  Precision 


Separation  in  Units  of  X/D 

Fig, 3 A comparison  of  results  obtained  by  current  an- 
alysis to  results  presented  in  [1]  and  [21,  Notice 
that  the  results  obtained  by  currant  analysis  are  much 
closer  to  those  obtained  in  [2],  where  an  approximated 
pulse  shape  was  used.  The  relative  measurement  preci- 
sion is  the  reciprocal  of  the  degradation  factor  pre- 
sented in  [21-  [4]  . 
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| Target  Separation  In  Units  of  X/D 

I • 

| 

f Fig. 5 A comparison  of  values  of  F^,  the  relative  separ- 

| • ation  measurement  precision,  obtained  by  current  analysis 

I to  those  presented  in  [11. 


tho  value  of  FA,  t’T  and  B‘t  can  differ  by  ag  much  as  3 orders 
of  magnitude!  This  severe  discrepancy  can  change  the  conslusiona 
regarding  the  capability  of  an  optical  sensor  system  drastically. 

in  Figures  7 through  9,  values  of  F^,  F^  and  Ffc  are  plotted 
against  target  separation  for  detector  widths  equal  to  0,31/D,' 

l.. 

2.44A/D  and  101/D,  respectively.  Although  the  value#  of  FA, 
and  Ffc  are  very  different  from  those  in  [1]  for  smaller  target 
separation,  the  target  separation  eorespondlng  to  F^-0.5  (3dB  de- 
gradation) does  not  change  a lot  due  to  the  numerical  errors. 
Figure  10  shows  the  value  of  target  separation  for  FA«*0.5  and  the 
"resolution  scale"  defined  in  [1)  as  a function  of  the  detector 
width.  It  indicates  that  regardless  of  the  numerical  problem 
presented  in  111  for  small  target  separation,  the  central  result 
of  111  concerning  about  the  oroposed  resolution  scale  is  still 
valid. 


IV.  SUMMARY  AND  CONCLUSION 

The  degradation  of  estimates  of  relative  target  strength, 
target  separation,  and  midpoint  position  for  a pair  of  closely 
spaced  objects  can  be  assessed  via  the  method  presented  here.  Wo 
also  show  that  the  results  oresented  in  11]  and  12],  which  have 
numerical  disagreements,  are  in  fact  theoretically  identical. 

The  cause  of  the  discrepancy  reported  in  ( 2]  is  due  to  numerical 
problems  occurring  in  a computer  program  used  by  Fried.  The  rms 
precision  with  which  we  can  measure  the  location  of  the  midooint 
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Fig. 7 Degradation  of  measurement  precision  on  the  target 
relative  intensity  as  a function  of  target  separation. 
Reeults  are  shown  for  detector  widths  equal  to  0.31/D, 
2.441/D  and  101/D. 


and  the  separation  of  the  pair  of  objects,  as  wdll  as  their  relative 
intensity  are  in  fact  equal  to  the  associated  Cramer-Rao  bounds 
oresented  here.  ,, 

In  examining  the  numerical  results  presented  here,  it  should 
be  emphasized  that  the  "optimistic"  conclusion  drawn  by  Fried  is  no 
longer  true  for  small  target  separations . in  same  casetij  when 
the  aeparation  is  significantly  less  than  A/D,  we  find  the  numerical 
result  presented  in  111  can  differ  from  ours  as  much  as  3 orders 
of  magnitude!  This  large  deviation  in  the  prediction  of  rms 
precision  can  change  the  conclusion  regarding  the  capability  of  a 
■particular  sensor  syatem  drastically! 

For  separation  greater  than  VD,  the  above-mentioned  numerical 
problem  seeniB  have  a very  little  effect  on  the  numerical  results 
for  the  prediction  of  rms  precision-  Since  target  separations 
corresponding  to  1^=0. 5 are  in  the  region  slightly  less  than 
to  much  greater  than  A/d,  Fried’s  major  conclusion  of  |1]  about 
the  proposed  resolution  scale  is  still  valid. 

Although  the  results  presented  here  are  associated  with  a 
particular  pulse  shape,  namely,  a one-dimensional  approximation  of 
a two  dimensional  unobstructed-circular-anerture,  diffraction-limited 
optics  with  a sharply  delineated  rectangular  detector,  the  analytical 
results  can  be  easily  extended  to  any  particular  optical  system, 
providing  that  the  pulse  shape  of  the  output  to  a unit  strength 
point  source  is  available.  The  noise  model  considered  here  is  an 
additive  white  Gaussian  noise.  For  a more  complicated  senBor  noise 
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APPENDIX!  CRAMER- RAO  BOUND  AND  ERROR  ANALYSIS  GIVEN  IN  IT 


In  [1],  the  following  error  covariance  matrix  is  considered t 


V - cov[{t0,  6T,  5 A 1 


which  is  slightly  different  from  the  error  covariance  matrix  we 
considered  in.  this  report.  The  error  covariance  matrix  associated 
with  the  Fisher  information  matrix,  F,  isi 


where 


V*  » cov[6R,  6T,  fit  1 > F "JL, 


n.'  v\ 


R - (1-ljA) /(!+»» A)  . 


It  is  easy  to  show  that 


whore 


* T 

V * AV  A 


-u+ha>*  o o 


" :7; 


Therefore,  the  Fisher  information  matrix  associated  with  the  error 
covariance  matrix  V becomes 


J - A"T  F A"1.  <A6) 

Substituting  Bgs.  (11a) -(Ilf)  and  (15) -(17)  into  (A6),  the  elements 
of  J becomes: 

J11  “ ^33  " 2a2(8ir2){[G2(0)+G2(T)  ] + 1*A2  (Gj  (0)  -G2  (T)  ] }/NQ,  (A7a) 

J12  “ F23  " ”2a2(4n2){AG2(0) }/Nq,  (A7b) 

J13  - -(1+ljA)"2^  - -2a2(2ir)  (H1(T)  )/N0,  (A7c) 

J22  * F22  " 2a2(2ii2){(G2(0)-G2(T)  ] + >»A2  [Ga  (0 ) +G2  (T)  ] (A7d) 

-2(A*A0)  ( (l-l|A2]H1(Ti  >2}/N0 

J23  " - ( 1+^A > ~ 2 F12  - 2a2  (2fr)  (A*AQ)  AH1(T)/Nq,  (A7e) 

, .4  , * (A7f) 

J33  “ <1+*»A)  *7 ii  - 2aNA0A0)  [1-Gq(T)1  [1+Gg (T ) 1 /NQ.  (A7f) 

where  we  have  used  the  following  identity: 

Fo_1  - l(l+l*A2)  + (l-%A2)r  (T)  ]"1/<4M0)  - (A*A0)/(2Nq)  (A8) 
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and  < A_A_ ) is  defined  in  Eq.  (13)  of  [11.  Notioe  that  the  J matrix 
o o 

is  differed  from  the  S matrix  given  in  the  Appendix  of  [1]  (Eqs. 

(A. 17)  - (A.  22))  by  a constant,  namely 


J.  . - - rrr~  s.  , 

H hqa0  ij 


for  all  i,j. 


If  we  compare  the  (i,j)^  element  of  M matrix  of  (1]  (Eqs.  (43a)- 
( 4 3 f ) ) with  S^.,  we  find  that 


Mij  " “Vo  Sij 


for  all  i,j. 


(A10) 


Rewriting  Eqs.  ( 39a) - ( 39c)  of  [1]  in  terms  of  matrix  notation,  we 
aotually  have  the  following  expression 


V “ DMD 


(All) 


where  D is  the  inversion  of  S.  Apnlying  the  relationship  we  find 
in  (A10) , (All)  becomes 


V « -N  A*  S”1. 
o o 


Substituting  (A9)  into  (A12) , we  finally  have 


? * 

2a*A_  . 

v 2 j"1. 


(A13 ) 


) :;3  •j'V: 


£i£.r* 


We  have  shown  the  covarianoe  matrix  V obtained  in  t 1 1 for  the 
error  analysis  equals  to  the  aasoeiated  Cramer-Rao  bound,  J_1,  by 
a constant  2a4A0/A0,  It  becomes  obvious  that  the  degradation 
factor  (a  ratio  of  two  CRB's)  F^,  F^,  and  derived  from  |1]  and 
current  analysis  should  be  identical. 
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